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ABSTRACT 

We present upper limits on the 850 /im and 450 /im fluxes of the warm hyperlumi- 
\Q ■ nous (bolometric luminosity L hol > 10 13 L ) galaxies IRAS P09104+4109 (z = 0.442) 

C\) ; and IRAS F15307+3252 (z = 0.926), derived from measurements using the SCUBA 

bolometer array on the James Clerk Maxwell Telescope. Hot luminous infrared 
C^Q ■ sources like these are thought to differ from more normal cold ultraluminous infrared 

(Xboi > 10 12 Lq) galaxies in that they derive most of their bolometric luminosities 
from dusty AGNs as opposed to starbursts. Such hot, dusty AGNs at high redshift 
are thought to be responsible for much of the mass accretion of the Universe that 
is in turn responsible for the formation of the supermassive black holes seen in the 
centres of local galaxies. The galaxy IRAS P09104+4109 is also unusual in that it is 
a cD galaxy in the center of a substantial cooling-flow cluster, not an isolated inter- 
acting galaxy like most ultraluminous infrared galaxies. Previously it was known to 
have large amounts of hot (T > 50 K) dust from IRAS observations. We now show 
that the contribution of cold dust to the bolometric luminosity is less than 3 per cent. 
Most ultraluminous infrared galaxies possess large amounts of cold dust, and it is now 
known that some cooling flow cluster cD galaxies do as well. Yet this object, which is 
| an extreme example of both, does not have enough cold gas to contribute significantly 

to the bolometric luminosity. We outline physical reasons why this could have hap- 
pened. We then provide a discussion of stategies for finding hot dusty AGNs, given 
the limitations on submillimctre surveys implied by this work. 
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1 INTRODUCTION 



Two extreme ultraluminous infrared galaxies (ULIGs) at 
redshifts z < 1 that are thought to derive most of 
their bolometric luminosity from a dust-enshrouded active 
galactic nuclei (AGNs) are IRAS P09104+4109 and IRAS 
F15307+3252. These objects differ from optical quasars, 
even those with "infrared bumps" (e.g. Sanders et al. 1989) , 
in that they are very faint at rest-frame optical wavelengths; 
their B-band luminosities are at least an order of magni- 
tude lower than their 25-pim luminosities, for example. These 
kinds of objects have attracted much recent attention be- 
cause of the possibility that much of the mass accretion in 
the Universe responsible for producing the high local density 
of supermassive black holes (Magorrian et al. 1998) might 
have occurred in high redshift analogues of these two galax- 



ies (Haehnelt, Natarajan & Rees 1998; Fabian & Iwasawa 
1999; Salucci et al. 1999). 

The galaxy IRAS P09104+4109 was identified by Klein- 
mann et al. (1988) to be an ultraluminous IRAS galaxy at 
z — 0.442. It was found in a follow-up survey of IRAS 
60-/im sources that had non-stellar colours and no opti- 
cal counterparts on the Palomar Sky Survey. Its bolomet- 
ric luminosity is 1.4 x 10 13 Lq, which is extremely high, 
even for a ULIG (it would be classed as a hyperluminous 
infrared galaxy by Sanders & Mirabel 1996). Radio, opti- 
cal, near- and mid-infrared imaging, spectroscopy, and spec- 
tropolarimetry (Kleinmann et al. 1988, Hines & Wills 1993, 
Soifer et al. 1996, Green & Rowan-Robinson 1996; Taniguchi 
et al. 1997, Evans et al. 1998, Hines et al 1999) suggest that 
this galaxy is powered by an obscured AGN. A hard X-ray 
detection of IRAS P09104+4109 with BeppoSax (Frances- 
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chini et al. 2000) confirms that this object contains a pow- 
erful AGN. 

The spectral energy distribution (SED) of IRAS 
P09104+4109 is unusual in that it peaks between 10 and 
30 /jm, implying that the dust responsible for reradiating 
the energy from the obscured AGN must be very hot, in 
excess of 100 K. This is very different from what is seen 
in other powerful (although lower in luminosity than IRAS 
P09104+4109) dust-enshrouded AGNs like the "warm" 
ULIGs Markarian 231 and IRAS 08572+3915 (Sanders et 
al. 1988a). These have dust temperatures closer to 50 K 
and SEDs which peak longward of 50 /xm. The dust in 
IRAS P09104+4109 is also considerably hotter than the 
dust in many infrared QSOs, like Markarian 1014 (Sanders 
et al. 1988b). That the dust in IRAS P09104+4109 is so 
much hotter than the dust in other ULIGs hosting pow- 
erful AGNs is perhaps suggestive that the dust in IRAS 
P09104+4109 is being heated directly by the AGN whereas 
in the other ULIGs, much of the far-infrared luminosity 
comes from colder, more diffuse, dust heated by star forma- 
tion (Rowan-Robinson & Crawford 1989, Rowan-Robinson 
& Efstathiou 1993, Rowan- Robinson 2000). 

Additionally, IRAS P09104+4109 is a cD galaxy in a 
rich X-ray cluster that is undergoing a very substantial cool- 
ing flow, about 1000 M yr _1 (Fabian & Crawford 1995, 
Crawford & Vanderriest 1996). In this way it is very differ- 
ent from other ULIGs, which are normally isolated peculiar 
galaxies, often showing signs of a recent interaction (Sanders 
et al. 1988a, Clements et al. 1996, Murphy et al. 1996). 

So IRAS P09104+4109 is a very unusual ULIG on three 
grounds: (i) its bolometric luminosity is very high; (ii) its 
SED is peaked at short wavelengths, implying a much higher 
dust temperature; and (iii) it is not a peculiar galaxy, but 
a cD galaxy in the center of a massive X-ray cluster that is 
undergoing a cooling flow. 

It is interesting to speculate as to whether these prop- 
erties are related; in particular one might ask whether or 
not the cooling flow can fuel the AGN and in doing so gen- 
erate the high bolometric luminosity and dust temperature. 
A related question is: does IRAS P09104+4109 have any 
cold (T < 50 K) dust at all, or is all the dust in this galaxy 
extremely hot and directly associated with the AGN. Al- 
most all other ULIGs (Sanders & Mirabel 1996) and some 
high-redshift quasars (McMahon et al. 1999) possess cold 
gas. Some cooling flow clusters do too, as inferred from their 
SCUBA submillimetre detections (Edge et al. 1999) and CO 
line measurements (e.g. Bridges & Irwin 1998). Cooling flow 
clusters also tend to have star-formation in progress at their 
centers (Cardiel et al. 1998, Crawford et al. 1999), which 
may be fuelled by this cold gas. This star formation could 
could directly produce the dust seen by Edge and collabora- 
tors. In IRAS P09104+4109, we have the additional evidence 
that there exists a high column density of X-ray absorbing 
gas (Iwasawa, Fabian & Ettori 2001) which is likely to be 
cold. 

The galaxy IRAS F15307+3252 (z = 0.926; Cutri 
et al. 1994) is another hyperluminous infrared galaxy 
thought to be powered by a luminous dusty AGN (Hines 
et al. 1995). It has somewhat cooler IRAS colours than 
IRAS P09104+4109 (but still warmer than the majority of 
ULIGs), suggesting its dust is at a lower temperature than 
that of IRAS P09104+4109. This object also differs from 



IRAS P09104+4109 in that it does not exist in an X-ray 
cluster. In fact, this object has yet to be detected in X-rays 
(Fabian et al. 1996, Ogasaka et al. 1997). It is presumably 
not associated with a substantial cooling flow. Unlike IRAS 
P09104+4109, there is no a priori reason to expect cold gas 
in this galaxy. 

At high redshifts (z > 2), objects which contain dusty 
AGNs are certainly known, such as IRAS F10214+4724 
(Rowan-Robinson et al. 1991) and APM 08279+5255 (Irwin 
et al. 1998). However the discovery of both of these objects 
has been serindipitous, often helped by strong gravitational 
lensing (see Eisenhardt et al. 1996 and Lacy, Rawlings & 
Serjeant 1998 for IRAS F10214+4724 and Egami et al. 2000 
and references therein for APM 08279+5255). More system- 
atic searches are required if we are to address the issues 
highlighted in the opening paragraph of this section. Sur- 
veys done at X-ray wavelengths offer one approach, and have 
begun to uncover some hot dusty objects (Wilman, Fabian 
& Gandhi 2000). But these will not find X-ray-quiet dusty 
AGNs like IRAS F15307+3252. Submillimetre searches with 
current or future bolometer arrays offer another method, but 
these will not find dusty AGNs if all the dust is hot; unless 
there is some cold dust present, sources would then fall be- 
low the submillimetre confusion limit (Blain et al. 1998). 
The presence or absence of cold dust in such objects there- 
fore has important implications for cosmological surveys. 

We did not detect this cold gas in either object, de- 
spite taking very deep integrations. In Section 2, we de- 
scribe the observations. In Section 3, we present the re- 
sults, compute SEDs, compare IRAS P09104+4109 and 
IRAS F15307+3252 to other infrared-luminous AGNs, dis- 
cuss the physical processes that may be responsible for gen- 
crating a lack of cold gas in these objects, and outline the 
consequences for systematic surveys of high-redshift dust- 
enshrouded AGNs. In Section 4 we summarize. Throughout 
this work, we assume the following cosmological parameters: 
H = 65 kms" 1 Mpc-\ fitter = 0.3, Ov = 0.7. 



2 OBSERVATIONS AND DATA REDUCTION 

Both sources were observed with the Submillimetre Com- 
mon User Bolometer Array (SCUBA; Holland et al. 1999) 
on the James Clerk Maxwell Telescope (JCMT) on Mauna 
Kea. This instrument consists of two bolometer arrays si- 
multaneously viewing the same region of the sky through a 
dichroic beamsplitter, arranged in a hexagonal grid with ap- 
proximately one telescope-beamwidth separating the detec- 
tors. We used standard filters: the "short-wavelength array" 
corresponded to an effective wavelength of observation of Ao 
= 450 fim and the "long- wavelength array" to Ao = 850 /im. 

All observations were conducted in photometry mode, 
in which the source is positioned onto the central pixel of 
each array. The data takes the form of a time- dependent 
voltage reading from each bolometer while the telescope sec- 
ondary mirror chops the source in and out of the beam, with 
a throw of 60 arcsec. In photometry mode, one "integration" 
is defined as a cycle of 9 on-off chopped-pairs positionally off- 
set by 2 arcsec, for a total of 18 seconds of observation. Pri- 
mary and secondary flux calibration sources were observerd 
to convert the bolometer voltages into fluxes. The details of 
the observations are summarized in Table 1, where we list 
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Table 1. Summary of SCUBA Observations 



Source 


uate 


Number of 


lotai time 


< T850 nm > 


< T450 )im > 


r 850 fim 


r 450 Mm 






integrations 


on source 






















mjy 


mjy 


IRAS P09104+4109 


20 Jan. 1999* 


50 


900 s 


0.25 


1.3 


4.24 ± 3.36 


-19.4 ± 34.2 




18 Mar 2000** 


120 


2160 s 


0.26 


1.4 


1.44 ± 3.61 


-13.35 ± 41.11 




4 Apr. 2000** 


280 


5040 s 


0.39 


2.2 


3.87 ± 4.18 


-29.29 ± 62.79 


IRAS F15307+3252 


6 Feb. 2000** 


175 


3150 s 


0.45 


2.6 


3.21 ± 4.36 


4.55 ± 83.24 




4 Apr. 2000** 


70 


1260 s 


0.36 


2.0 


1.79 ± 4.13 


6.98 ± 36.23 



* Observations performed by K. Hodapp (University of Hawaii) 
** Observations in service mode. 



integration times, the local atmospheric optical depth at the 
time of observation, and the measured fluxes. 

The alternative positive-and-negative waveform of the 
chopping secondary was taken out of the voltage stream. 
Next, the time-dependent variations of the sky brightness 
were estimated using the signal variations in the bolome- 
ters at the periphery of each array, and then were sub- 
tracted from the data streams. Two methods of estimating 
the sky variations (based on using either the mean or the me- 
dian of the outer bolometer readings) were employed, which 
proved mutually consistent to within 15% in all cases. We 
then subtracted the mean of these two sky brightness esti- 
mates from the values at the central (on-source) pixel. In all 
cases, this sky subtraction reduced the uncertainty on the 
mean source voltage while the improved value remained con- 
sistent with the non-sky-subtracted measurements within 
the latter's larger uncertainties. For each dataset (observ- 
ing run/source combination), the mean and the variance of 
the central pixel voltage over time was then determined. 

This entire procedure was first performed on the obser- 
vations of the secondary flux calibrators to determine the 
flux conversion factor (FCF), by which the mean voltage (in 
mV) is converted into the observed source flux (in mjy). 
The following sources were used as flux calibrators: Mars, 
CRL618, IRAS+16293, OH231.8, and IRC+10216. Their 
assumed fluxes were taken from the SCUBA calibration 
web page (http:www.jach.hawaii.edu/JACpublic/JCMT/ 
continuum-observing/SCUBA/). The FCFs were deter- 
mined independently for each array (i.e. wavelength) and for 
each observing run. Most of the data were acquired in ser- 
vice mode under "Category 3" weather conditions, which are 
considered suitable for observations using the Long Wave- 
length Array but inadequate for use with the Short Wave- 
length Array. Appropriately, the 850-/im FCFs were mea- 
sured to be more stable than the 450pm FCFs. During a 
single observing run, the 850/im FCFs typically varied by 
~10 — 15%, and showed run-to-run variations of ~30%. 
The 450-/im FCFs varied much more widely, up to a fac- 
tor of two or more, both during a single observing run, and 
between runs. The variation of the calibrator FCFs was one 
of the most significant contributions to the (systematic) un- 
certainty in our results. Faint-object photometry at submil- 
limetre wavelengths is a difficult and error-prone exercise, 
particularly under below-average weather conditions such as 
these, and the uncertainty in our final fluxes (see Table 1) 



is large. Neither IRAS P09104+4109 nor IRAS 15307+3252 
were detected within our errors: only upper limits to the 
source fluxes were obtained. 

In addition to these observations, an upper limit to 
the 350-/im observed-frame (243-/im rest-frame) flux of 
IRAS P09104+4109 was obtained at the Caltech Submil- 
limeter Observatory in October and December of 1997 
using the SHARC camera (Hunter, Benford & Serabyn 
1996). Although superior weather conditions were experi- 
enced, the collecting area, observing efficiency and stabil- 
ity of the SHARC camera are less than that of SCUBA. 
This is primarily because the SHARC detector array is one- 
dimensional and has several inoperative elements near the 
array center. The upper limit to the flux obtained (~300 
mjy) does not contribute significantly to the 450 and 850 
limits discussed here. 



3 RESULTS AND DISCUSSION 

Combining data from all the runs listed in Table 1, we find 
the following fluxes: for IRAS P09104+4109, 3.18±2.12 mjy 
(3(J < 9.54 mjy) at 850 ^m and -18.77 ± 24.25 mjy (3a < 
72.75 mjy) at 450 /im, and for IRAS F15307+3252, 2.46 ± 
3.00 mjy (3cr < 11.46 mjy) at 850 ^m and 6.59 ± 33.22 mjy 
(3a < 106.3 mjy) at 450 /jm. For IRAS F15307+3252 the 
constraints from our results are comparable to those from 
Yun & Scoville (1998), although at shorter wavelengths. 
SEDs are plotted in Figures 1 and 2. 



3.1 The lack of cold gas in IRAS P09104+4109 

It is clear from Figure 1 that IRAS P09104+4109 has very 
little cold dust; reradiated energy from a cold dust compo- 
nent contributes a fraction below 3.0 % of the bolometric 
luminosity. The contribution to the 850-/im flux that we 
would expect from the hot dust responsible for the 25-pim 
IRAS detection accounts for only 1 /400 of our observed 2a 
upper limit; if this was the only dust in the galaxy, it would 
take several weeks to detect with SCUBA! 

The absence of cold dust in a powerful AGN like this is 
certainly unusual. High-redshift optical quasars often have 
cold dust (McMahon et al. 1999), so one might expect an in- 
frared quasar like IRAS P09104+4109 to have some too, par- 
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8 10 12 14 16 18 20 

Log 10 (iVHz) 

Figure 1. The spectral energy distribuion of IRAS 
P09104+4109. The points arc: open circles - 0.44 /im, 0.55 jim, 
0.64 (im, 0.79 fim, 3.8 fim, 10.1 fim and 20 (all in the ob- 
server frame) from Klcinmann ct al. (1988); open squares - radio 
continuum at 6 cm and 20 cm from Hincs & Wills (1993); filled 
circles - 1.25 fim, 1.65 fim and 2.2 fim from Soifer et al. (1996); 
open triangles - 12 fim, 25 lira, 60 (im and 100 fim limit from 
the IRAS Faint Source Catalog; filled triangles - 450 fim and 850 
fim SCUBA 2a limits from this work; filled squares — 0.5-7 keV 
Chandra and 14—40 keV BeppoSax detections from Iwasawa et 
al. (2001). The solid line is the best-fitting modified (by a dust 
cmissivity k oc i/ 15 law) blackbody fit to the IRAS detections; 
this has a dust temperature T = 206 K and peaks at 13 fim, but 
is a statistically poor fit to the three data points. The solid line 
is the maximum contribution from a cold (37 K) dust component 
given the limits imposed by our SCUBA data. 

ticularly since the X-ray measurements (Iwasawa et al. 2001) 
tell us that plenty of cold gas is present. 

The solution to this apparent paradox is presumably 
tied up in the way that gas condenses out of the hot X-ray 
phase in a cooling flow. If gas condenses in such a way that 
there is no new star formation, dust generation will be inef- 
ficient - carbon atoms associated with this cooling gas are 
more likely to interact with four hydrogen atoms and become 
methane rather than with numerous other carbon atoms and 
become dust grains. This could explain the presence of large 
amounts of cold X-ray absorbing gas (Iwasawa et al. 2001) 
and no cold dust. If the cooling flow is also responsible for 
fuelling the AGN, any star formation that happens along 
the way must be confined to regions very close to the AGN. 
Red giants and supernovae associated with this star forma- 
tion could generate the hot dust seen by IRAS. 

3.2 The lack of cold gas in IRAS F15307+3252 

We also did not detect cold dust in IRAS 15307+3252. In 
this case the upper limit on the fraction of the bolometric 
luminosity that could come from a cold dust component is 
0.3 %. 
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Figure 2. The spectral energy distribuion of IRAS F15307+3252 
The points are: open circles - 0.44 fim and 0.70 fim, (all in the 
observer frame) from Cutri et al. (1994); open square - radio 
continuum at 21 cm from the FIRST survey of Becker, White & 
Helfand (1995); filled circles - 1.25 fim, 1.65 fim and 2.2 fim from 
Soifer et al. (1994); open triangles - 12 fim, 25 fim, 60 fim and 
100 fim limit from the IRAS Faint Source Catalog; filled triangles 

- 450 fim and 850 fim SCUBA 2a limits from this work; filled 
pentagon — 650 fim lmit from Yun & Scovillc (1998); filled squares 

- 2 -10 keV ASCA limit from Ogasaka ct al.(1997). The solid line 
is the best-fitting modified (by a dust emissivity k oc i/ 1 - 5 law) 
blackbody fit to the IRAS detections; this has a dust temperature 
T = 126 K and peaks at 21 fim. The solid line is the maximum 
contribution from a cold (37 K) dust component given the limits 
imposed by the data of Yun & Scovillc (1998). 

While we had no a priori reason to expect dust in thie 
galaxy, this is still a suprising result. High-redshift optical 
quasars often have cold dust (McMahon et al. 1999). Yet 
this galaxy, which is opically faint (presumably due to in- 
ternal extinction) and is such a strong 60- fim emitter, has 
very little. This result is in some ways more perplexing than 
the non-detection for IRAS P09104+4109 since the IRAS 
colours are cooler. 

These two non-detections have implications for the dis- 
covery of hot, dusty AGNs at high-redshift using infrared 
and submillimetre surveys. We investiate these implications 
in the rest of this section. 

3.3 Application to high redshift 

In Figure 3, we plot an optical-infrared vs. infrared- 
submillimetre colour-colour diagram for the sample of dust- 
enshrouded AGNs presented in Table 1. In the rest of this 
section we will use this diagram to assess the efficacy of find- 
ing dusty AGNs at high redshift by a variety of methods. 
Some particular features worth noting are: 
(i) higher temperature sources are to be found near the top 
of this diagram, but this does not tell us immediately that 
all sources in this region of the diagram are AGN-powered. 
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For example, IRAS F10214+4724 is to be found near the 
top of the diagram yet it may derive a substantial fraction 
of its total luminosity from an obscured starburst (Green & 
Rowan- Robinson 1996); 

(ii) how far to the right of this diagram is a measure of 
how optically bright an object is, which in turn has impli- 
cations for our ability to infer the presence of a powerful 
AGN using optical spectroscopy. For example, the "Ql", 
"Q2", "CLOV", and "APM" points all represent objects 
with AGNs so powerful that they likely are responsible for 
the high bolometric luminosities of the objects; 

(iii) this diagram gives no indication of the absolute lumi- 
nosities of the objects, rather the relative luminosities at dif- 
ferent wavelengths. This has implications for the detectabil- 
ity of the objects. For example, APM 08279+5255 is near the 
top of this diagram yet has been detected at submillimetre 
wavelengths with SCUBA (Lewis et al. 1998) ; this is becuase 
both its bolometric luminosity and gravitational magnifica- 
tion are so high that its submillimetre flux is still measur- 
able. Conversely, IRAS P09104+4109, as we saw here, was 
not detected, although its spectral energy distribution is per- 
haps not all that different from that of APM 08279+5255 
(note that the increased distance of APM 08279+5255 rela- 
tive to IRAS P09104+4109 is compensated for by the nega- 
tive K-correction (see Section 3.3.1); 

(iv) Radio and X-ray luminosities can be quite different for 
objects in similar regions on Fig. 3. This can also have impli- 
cations regarding the detectability of the different objects. 

3.3.1 Submillimetre surveys 

The main result of the present paper is that the hot dusty 
AGNs IRAS P09104+4109 and IRAS F15307+3252 have 
very little if any cold gas. So if all of the dust-enshrouded 
accretion in the Universe took place at high redshift in 
objects like these, we would not see them in submillime- 
tre surveys like those carried out with SCUBA (e.g. Smail, 
Ivison & Blain 1997, Hughes et al. 1998, Cowie & Barger 
2000). The only ones we would see in such surveys are 
(1) extremely luminous, possibly highly-magnified, objects 
like APM 08279+5255, and (2) the most distant objects 
at very high redshifts z > 6: for sources with dust emis- 
sivity k v ~ v 1 ' 5 (e.g. Blain et al. 1999), the observed flux 
~ (1 + z) 4 " 5 <Il(z)~ 2 where cIl(z) is the luminosity distance, 
so that the negative K-correction more than compensates 
for the dimming due to the increased distance by a signifi- 
cant amount (see Fig. 4). However, we do not expect the vast 
majority of sources to fall into either of these two categories. 

3.3.2 Infrared surveys 

Mid- and far-infrared survey s, like those to be carried ou t 
with the SI RTF satellite ( frittp: / /ipac. caltech.edu/sirtf/ ) 
offer a far more productive strategy for finding hot dusty 
AGNs at high redshift: both IRAS P09104+4109 and IRAS 
F15307+3252 would be easily detected at redshifts z < 4. 
The absence of submillimetre detections at the positions of 
bright mid-infrared sources would suggest the presence of 
hot dusty objects. Establishing that such objects are AGN- 
powered will be more difficult, particularly if they are op- 
tically too faint (as both IRAS P09104+4109 and IRAS 



F15307+3252 would be at high z) to permit spectroscopy 
with the largest ground-based telescopes. However, the anal- 
ogy with the lower-redshift objects studied here in the con- 
text of the observations and models outlined in Section 1 
would be suggestive. 

3.3.3 Optical, radio and X-ray surveys 

Surveys at other wavelengths can also uncover hot dusty 
AGNs, but in each case the bias introduced by the selec- 
tion technique is severe and complete samples cannot be 
obtained. The reason for this is that so little of the bolomet- 
ric luminosity of hot, dusty sources emerges outside the far- 
and mid-infrared wavebands. 

Optical surveys can uncover extreme objects like APM 
08279+5255 (Irwin et al. 1998). But these represent the rare 
instances where the AGN is visible directly at rest-frame ul- 
traviolet wavelenghts, possibly because it has burnt through 
its dust shroud. There still must be enough dust present, 
however, to allow it to be red enough to be identified in 
these surveys. It must also be extremely luminous, and in the 
case of APM 08279+5255, highly magnified. The majority of 
dusty AGNs at high redshift may well have properties very 
substantially different from this and would consequently be 
missed in optical surveys. Both IRAS P09104+4109 and 
IRAS F15307+3252 would be missed in optical surveys, 
for example; their optical colours and magnitudes are un- 
remarkable. 

Radio surveys can also reveal some dusty AGNs like 3C 
318 (Willott, Rawlings & Jarvis 2000). But, as for optical 
AGNs (Goldschmidt et al. 1999), it seems likely that the 
majority of such objects are not radio-loud. 

Addtionally X-ray surveys would certainly uncover 
some hot, dusty AGNs like IRAS P09104+4109 (although 
not IRAS F15307+3252). The presence of a population of 
sources seen by Chandra but not SCUBA (Fabian et al. 2000; 
see also the discussion of energy budgets and the relation to 
the local MDO density of Trentham & Blain 2001 and the 
more recent observations of Bautz et al. 2000, Hornschemeier 
et al. 2000, and Barger et al. 2001) may have revealed such 
a population already. Infrared observations of these sources 
with SI RTF will tell us whether or not these objects contain 
substantial amounts of hot dust. If so, the analogy with the 
local objects IRAS P09104+4109 and IRAS F15307+3252 
outlined in Section 3.2.2 would suggest that they might be 
AGN-powered. Two X-ray sources have already had ISO 
mid-infrared counterparts identified (Wilman et al. 2000), 
but longer wavelengths are required to establish the dust 
temperatures. 



4 CONCLUSION 

In summary, we have not detected cold, dusty gas in ei- 
ther IRAS P09104+4109 or IRAS F15307+3252. Dust with 
a temperature of 40 K, which is what dominates the bolo- 
metric luminosity in most local ULIGs and in the high- 
redshift SCUBA sources recently discovered in various sur- 
veys (e.g. Blain et al. 1999) contributes at most 3 % of the 
bolometric luminosity of these objects. 

These results were somewhat surprising given that 
IRAS P09104+4109 is associated with a massive cooling flow 
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Figure 3. Optical (B) - infrared (25 fim) vs. infrared (25 (im) - submillimctrc (850 fim) colour-colour diagram for 
the objects listed in Table 1. All quantities refer to the rest-frame. The 850-/im luminosities are computed assuming 
a modified Raylcigh-Jcans spectrum with dust cmissivity k ~ i/ 1 - 8 (e.g. Blain et al. (1999)). The 25-/iin luminosities 
are computed from IRAS measurements, using logarithmic interpolation and/or extrapolation (but never beyond 
130 fim) if necessary. The B-band luminosities are computed from K-corrected optical and near-infrared fluxes. The 
symbols and labels refer to the objects listed in Table 2. Additional comments are: "P09" - The dotted line represents 
the permitted values given a lower limit to the 850 fim flux equal to the contribution from the hot blackbody that 
is inferred from the IRAS 60 fim and 25 ^m detections — the solid line in Fig. 1; "Q3", "Q4", "SI", "S2", "S3", 
"HR10", "4C" and "8C" - The dotted line represents the locus of permitted values given the absence of a mid- 
infrared rest-frame detection. The upper-limit to the 25— (im luminosity comes from IRAS measurements, the lower 
limit from an extrapolation of the submillimctrc SED to mid-infrared wavelengths assuming a modified Raylcigh- 
Jeans spectrum as above with dust temperatures of 40 K for submillimetre galaxies (Blain et al. 1999; Trentham et 
al. 1999) and 50 K for quasars (McMahon ct al. 1999), and the open circle (where present) is derived from an ISO 
15-fim detection, converted to a rest-frame 25— fim flux assuming a power-law spectrum with index a = —1.7 (Blain 
et al. 2000); "WFG" — the mid-infrared luminosity is similarly derived from the IRAS 15-fim detection, and the 
assumed rcdshift is that derived by Wilman et al. (2000) using the HYPERZ population synthesis code (Bolzoncllo, 
Pello & Mirallcs 2000). For "Ql", "Q2" and "Q3", median values are plotted. For the "S88a" and "S88b" sequences, 
each point (the open or filled triangles) represents a galaxy, ordered as in those papers (the large point represents 
the first galaxy i.e. Arp 220 for "S88a" and IRAS 12071-0444 for "S88b") and connected by the lines. For IRAS 
F10214+4724, APM 08279+5255 and the Cloverleaf quasar, the positions on this diagram have not been corrected 
for differential magnification (these are strongly lensed sources) and corrections will shift the points downwards and 
to the left if the sizes of the emitting regions increase with longer wavelength. For the weakly lensed soures (the 
three SCUBA galaxies), differential magnification (Blain 1999) is a small effect. 
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Table 2. Objects plotted in Figure 3 



Notation Object + Description 

(additional to text) 



Redshift(s) 



Data references 



P09 IRAS P09104+4109 

F15 IRAS F15307+3252 

F10 IRAS F10214+4724 

APM APM 08279+5255 

CLOV H1413+117 

"Clovcrleaf Quasar" 

Ql Low-z + high L PG quasars 

L bol > 10 12 L 

Q2 Low-2 + low L PG quasars 

L bol < 10 12 L Q 

Q3 High-z quasars 

Q4 SMM J04135+1027* 

Submillimetre-selected quasar 

51 SMM J02399-0136 

52 SMM J14011+0252 

53 SMM J02399-0134 

HR10 ERO J164502+4626.4 

(Object 10 in Hu & Ridgway 1994) 

3C 3C318 
4C 4C41.17 
8C 8C1435+63 

CF1 A2390 cD 

cooling flow cluster 

CF2 A1835 cD 

cooling flow cluster 

WFG CXOUJ215333. 2+174209 

(Object B in WFG00) 

S88a Ultraluminous infrared galaxies 

S88b Warm ultraluminous galaxies 



0.44 
0.93 
2.29 
3.87 
2.56 

0-2 (mostly < 0.5) 

0-2 (mostly < 0.5) 

3.9 - 4.7 
2.83 



2.80 
2.56 
1.06 

1.44 



1.57 
3.80 
4.25 

0.23 

0.25 

0.575 



< 0.1 

< 0.2 



Klcinmann ct al. 1988, IRAS Faint Source Catalog, this work 
Cutri et al. 1994, IRAS FSC, this work, Scovillc 1997 

Rowan-Robinson et al. 1993 
Irwin et al. 1998, Lewis et al. 1998, Ibata et al. 1999 
Barvainis ct al. 1995, Hughes ct al. 1997 

Sanders et al. (1989) 

Sanders et al. (1989) 

McMahon et al. (1999) 
Knudsen et al. (2000) 



Ivison et al. 1998 
Ivison ct al. 2000 
Soucail ct al. 1999, Smail ct al. 2000 

Dey et al. 1999 



Willott et al. 2000 
Graham ct al. 1994, Archibald et al. 2000 
Lacy ct al. 1994, Archibald et al. 2000 

Edge et al. 1999 

Edge et al. 1999 

Wilman et al. 2000 (WFG00), Lemonon ct al. 1998 



Sanders ct al. 1988a 
Sanders et al. 1988b 



in which there is independent evidence for cold gas (Iwa- 
sawa et al. 2001). Furthermore, cold dusty gas has been de- 
tected in a number of optically bright quasars (McMahon et 
al. 1999) so it was unexpected that it would not be detected 
in these two objects which are optically faint, presumably 
due to internal extinction. Some mass of cold dust may still 
be present, and may be responsible for some of this internal 
extinction, but it contributes negligibly to the bolometric 
luminosity. 

These two hot dusty galaxies are thought to derive most 
of their bolometric luminosity from AGNs as opposed to 
starbursts, based on a number of detailed observational re- 



sults and models. Dust-enshrouded AGNs that are not seen 
in optical surveys are thought to be responsible for much of 
the mass accretion in the Universe that produces the high 
local density of supermassive black holes seen in the cen- 
tres of nearby galaxies. Our results suggest that the best 
way to find such obejets is to look for powerful mid- and 
far-infrared emitters in SI RTF surveys that have very low 
submillimtre SCUBA fluxes. Some such objects may have al- 
ready been seen in joint Chandra-SCVBA surveys, but there 
could also exist a population of hot dusty AGNs (like IRAS 
F15307+3252) that are not X-ray-loud. 
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